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Abstract 
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new light. 
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1 Introduction 



One puzzle of the random matrix theory approach jl], g, |3| to the computation of Dirac operator 
spectra has been the case of one fermion species. While the (massless) microscopic spectral densities 
associated with the three chiral analogs of classical matrix ensembles have been found to be different 
|TJ , the corresponding finite- volume partition functions are all equal [Q, ||[]. This means that all three 
microscopic spectral densities, despite being very different, should lead to the same spectral sum rules 
of, for instance, the kind [y indicates the (positive) topological charge): 

in 




4(i/ + 1 



which indeed they do This holds also in the case of double- microscopic spectral densities, where 
fermion masses are kept finite (and scaled with volume at the same rate as the eigenvalues) ||. 

As we shall show in this paper, there are in fact simultaneously two reasons for why three different 
universality classes of gauge theories can share the same finite- volume partition functions, have in 
common an infinite set of spectral sum rules, and yet have very different microscopic spectral densities. 
Let us first focus on the chiral unitary ensemble, which corresponds to SU(N C ) gauge theories with 
N c >3 and Nf fermions in the fundamental representation of the gauge group. As was explained in 

ref. 0, to compute the (double-) microscopic spectral density pg (£; fix, . . . , UNf) corresponding to 
this random matrix ensemble, one needs, in addition to the finite- volume partition Z^ N f\fi\, ... ,p,j^ f ) 
of the theory with Nf flavors, also that of a theory with two additional fermion species (of imaginary 
mass), Z( N f +2 \/j,i, . . . , /ijv/j ^C; *C)- So for Nf = 1 we need also the partition function for Nf = 3 (and 
the finite-volume partition functions corresponding to the two other random matrix ensembles for 
Nf = 3 are different). In addition, and this is one of the main points of this paper, also the precise 
relationship between the (double-) microscopic spectral densities and the corresponding finite-volume 
partition functions turn out to be different. 

One case which stands on a rather special footing is that of the (non-chiral) unitary ensemble. It 
has been conjectured to be relevant for QCD-like theories in (2 + 1) dimensions with an even number 
of fermions species. In section 2 of this paper we begin by deriving the relations which allow us to 
directly compute all double-microscopic spectral correlators from the finite- volume partition functions 
of QCD3 alone. We also show how the universal limits of the orthogonal polynomials of the matrix 
model formulation can be computed directly from a QCD3 partition function. In section 3 we turn to 
the chiral unitary ensemble described above, where we do the extension of the analysis of ref. to the 
case of arbitrary topological charge u, and also here demonstrate how the universal microscopic limit 
of the associated orthogonal polynomials are given in terms of finite-volume partition functions. In 
section 4 we describe the extent to which we have been able to derive analogous relations for the two 
remaining categories of (chiral) random matrices: the orthogonal and symplectic ensembles. We end 
in section 5 with our conclusions and a discussion of the notion of matrix-model universality || ^, ^] . 
as seen in this new light. 



1 When considered in the appropriate "mesoscopic" scaling region. Finite-volume partition functions are always taken 
to be in this regime in what follows. 
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2 The Unitary Ensemble 



The unitary ensemble has been conjectured to describe, in the microscopic limit, the Dirac operator 
spectrum of (2 + l)-dimensional SU(N C ) gauge theories with N c >3 and an even number of fermions 
Nf Q. The relevant partition function of that ensemble is, for massive fermions, 

N f 

Z( N f\m u ...,m Nf ) = fdM J] det (M + im f ) e ~ NtTV ^ , (2) 

where the integration is over the Haar measure of hermitian N x N matrices M. The even potential 
V(M 2 ) has been left unspecified, and we write it in general as 

V(M 2 ) = ]T %M 2k . (3) 
fc>i 2k 



It is important to note that in order for this to possibly describe the microscopic limit of the above 
class of gauge theories, the even number of fermions have been regrouped into two sets of which one 
half is assigned the same values as the other half, except for a minus sign, i.e., the mass matrix can 
be taken to be of the form 

diag(mi,m 2 , . . .,m Nf / 2 , -mi, -m 2 , • • • , -m Nf/2 ) ■ 

With this assignment, the fermions can effectively be regrouped into 4-spinors, the dynamics of which is 
reminiscent of the corresponding (3+l)-dimensional gauge theories. The associated "chiral symmetry", 
which really is a flavor symmetry, has been conjectured to break spontaneously according to U(NA — > 
U(Nf/2) x U(Nf /2) (recall that Nf is even) 0|. An order parameter for this symmetry breaking 
pattern is the absolute value of the chiral condensate £ =Ya I witfill/Nf- I n matrix model terminology 
a non-zero condensate translates into a non-zero spectral density at the origin, p(0) ^ 0. For the 
underlying field theory, this is simply a (2+l)-dimensional generalization of the Banks-Casher relation 
between the chiral condensate and the spectral density of the Dirac operator, evaluated at the origin. 

In terms of the eigenvalues Aj of the hermitian matrix M we have, ignoring all irrelevant overall factors, 

roo N / N f /2 \ 

m\ mi ,...,m Nf ) = / n I** n( A ?+™/) e ~ w(Af) • w 

J -°°i=i V f=i J 

The imaginary part arising from the term inif in eq. (§) has disappeared due to the pairwise grouping 
of masses discussed above, and the final expression is, except for the "fermion determinant", a con- 
ventional eigenvalue integration in the unitary ensemble with arbitrary even potential. The massive 
spectral correlators of this ensemble have recently been derived in the double-microscopic limit where 
masses and eigenvalues are considered on the same scale of magnification around the origin M, and 
have been proved to be universal in the random matrix theory context . 



We are now in a position to make use of a nice result from ref. ]TT| (see also ref. [12]), where the 
spectral two-point function, the kernel Kn(x, y), was expressed in terms of a slightly modified random 
matrix integral. By definition, 

N f /2 at_i 

K i ^'' / \x,x';m 1 ,...,m Nf ) = e"T (V(* 2 )+V(x»)) jj ^( x 2 + m 2 )(x / 2 + m 2) £ P i {x)P i (x') , (5) 

/ »=0 
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where P n (x) is an nth order polynomial orthogonal with respect to the even weight function 

-NV{x 2 ) 



WiX) 



Nf/2 

n(x 2 + 



m 2 ) e 



(6) 



and this polynomial is hence also a function of all masses rrii. For convenience we consider here the 
normalization where the polynomials are chosen orthonormal on the real line. By making use of the 
orthonomality one can write the kernel in the form of a matrix integral: 



K N (x,x ;mi, . . .,m Nf ) 



e 2 



K(y{x 2 )+V{x l2 )) u N f /2 



Uf J(x 2 + mj)(x' 2 +m 2 f ) 



Z( N f){m\i . . .,m Nf ) 

N f /2 



roc N-lf N f /2 

/ Y[\dXi(Xi-xX\i-x')l[(X 2 + m 2 f )e 

J -°° i=i \ /=i 



-W(A?) 



detijXy- 1 



'•(7) 



Although the last matrix integral is over (JV — 1) eigenvalues only, this distinction becomes irrelevant 
in the large-iV limit, which we will consider below. This gives us a very convenient expression for the 
spectral two-point function in the large- A r limit: 



K N (x,x ;mi, . . .,m Nf/ 



e 2 



(V(x 2 )+V{x' 2 )) 



Nf/2 



II J{x 2 + m 2 )(x' 2 + m 2 ) 



,m Nf ,ix,ix') 
Z( Ar /)(mi, . . .,m Nf ) 



(8) 



The partition functions are of course symmetric in the mass entries, and the position of the additional 
fermion masses in the argument list is therefore immaterial. 

In writing the kernel in the form (|8|) we have made use of the observation that the insertion of the 
additional factor of (Aj— x)(Aj— x') in the matrix integral can be viewed as considering a theory with two 
additional fermion species, of imaginary masses (cf. eq. @). Except for the shown prefactors, we have 
thus expressed the two-point spectral correlation function in terms of the ratio of two unitary-ensemble 
matrix model partition functions. All higher-order correlation functions are then also manifestly 
expressed in terms of these partition function through the usual factorization relation of the large- N 
limit. Finally, also the conventional (macroscopic) spectral density is expressed in such a manner, now 
in terms of two additional fermion species degenerate in (imaginary) masses: 



x; mi, . . . , mjv 



lim K^ f (x, x; mi, . . . , mjv f ) 

N— >oo 



(9) 



The above results hold in all generality in the planar limit. As in ref. ]7|, we now consider specifically 
the double-microscopic limit in which 

C = Trp{0)Nx and = 7rp(0)Nmi (10) 

are kept fixed in the limit N — >oo. If we identify S = 7rp(0), this is precisely the mesoscopic scaling 
region of the finite volume partition functions of the underlying gauge theories. In ref. § these 
partition functions were argued to be describable in terms of a very simple chiral Lagrangian, in 
analogy with the Leutwyler-Smilga analysis in (3+1) dimensions: 

^(/i!,...,^) = J dU exp[Tr(MUT 5 U^} . (11) 
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Here A4 is the mass matrix discussed above, rescaled by NT,: 

M = dia,g((ii,H2,...,fi Nf /2,-fii,-H2,...,-iiNf/2) > ( 12 ) 

and r5 = (ljy /2, — ljvy/2) where ljvy/2 is an (Nf /2)x(Nf /2) unit matrix. The above partition function 

can be explicitly evaluated by making use of the Harish-Chandra-Itsykson-Zuber integral [13]. The 
result is ||: 

/ A({m}) A({-m}) \ 

z^'Hm) = V A{M) ' > ( 13 ) 

where the (Nf/2)x(Nf/2) matrix A({yUj}) is defined by 

= (/if)^- 1 ^ . (14) 
The denominator is given by the Vandermonde determinant of rescaled masses: 



N 



f 

MM) = II^-Mi) • (15) 

i<j 

In the mesoscopic scaling limit the matrix model partition functions Z^ N f\^i\, . . . , /zjvy) should equal 
the field theory partition function 2>( N f'(j/,i, . . . , /ijV/) up to an irrelevant (mass-independent) normal- 
ization factor. Furthermore, in this scaling regime the prefactor of exp[— (N/2)(V(x 2 ) + V(x' 2 ))] in the 
expression for the kernel @ becomes replaced by unity. For the kernel this leads us to the following 
master formula: 



(16) 



k s (ccj/xi,...,^) = c n v/(C 2 + ^)(C ,2 + ^) — Wf^^T^t) • 

Similarly, the double-microscopic spectral density becomes 

and all other microscopic correlators are given by 

PS {.Cl, ■ ■ ■ Xn, fil> ■ ■ ■ , m f ) = de t K s(Ca,(b',Vl, ■ ■ ■ ,fJ>N f ) ■ ( 18 ) 

a,b 

In these expressions there is still one overall normalization factor, C, which remains to be fixed. As 
discussed in ref. J5J, the simplest way to fix it is to make use of the matching condition between the 
microscopic spectral density and macroscopic spectral density, 

lim pf f) ((;»!,..., fi Nf ) = p(0) = (19) 

where we have inserted the conventional normalization [l| (this conveniently makes \i{ = Nrrii and 
(i=N Zi ). 



To illustrate the power of the master formula fll6|), consider the simplest case of quenched fermions, 
which formally corresponds to Nf = 0. In this case the origin A = is not singled out, and the 
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kernel Ks(C, (') should reduce to the famous sine-kernel of the unitary ensemble. Indeed, the relevant 
partition function for two flavors in this case reads, from eq. fll3|), 

Mi - M2 

while the finite- volume partition function for Nf = is a trivial constant, which we set to unity. The 
kernel is therefore given by 

KsK.C) = C 2ai f_- ( \ (21) 
and the microscopic spectral density, as expected, becomes just a constant: 

P P(0 = 2C . (22) 



The matching condition (IS) hence gives C = l/(27r), and thus 
- a novel derivation of the sine-kernel. 

Because we know the analytical form of the finite-volume partition function for any number Nf of 
(massive) fermion species, we can immediately write down the general expressions for all double- 
microscopic spectral correlators. From eq. (13) it follows that 

/B({^,Ci}) B({-^,-Ci}A 

-j£ aet \B({-^,( 2 }) B(fo,-C 2 » ) 

s (Ci, G; Mi, • • • , m) (Ci _ + (c| + de /A( {M ) A( { - M )\ 

V A ({-^i) A ({M) J 

(24) 

where A({^j}) is the (Nf/2) x (Nf/2) matrix defined previously in eq. (14), and B({/Zj,£j}) is an 
(JV//2+l)x(JV//2+l) matrix defined by 

Nt N f 
B ki = A M for 1 < k < -f ; 1 < I < -f + 1 

Nf Nt 
Bki = (iCi,2)'-V^ for k = ^- + l ■ l<l<-± + l t (25) 

where in the last line the entry is either £i or (2, as indicated explicitly in eq. (| 
It similarly follows from eq. ( |i~7|) that the double-microscopic spectral density is 

/ B({ Mi ,C» B({- Mi ,-C}) 

wv , -c et {c({-^Q) c({^,-c» , 



aet {A({-^}) A({ Mi }) y 

where C({^,C}) and C({^,C}) are (N f /2+l) x (JV//2+1) matrices defined by 

iV> N f 
C kl = C kl = A kl for 1 < k < -j- ; 1 < I < -j- + 1 
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C ki = -C M = «)'- 2 e i? (iC + /-l) fOT ^ = ^ + 1 ; 1 < I < ^ + i . (27) 

All other double-microscopic spectral correlators are then also known explicitly, using the relation 
©■ 

It remains to fix the overall constant C in this general case. We again do it by the matching condition 
(U). This gives 

C= — . 

All double-microscopic spectral correlators, for any even value of Nf, are then completely determined. 
Recently these double-microscopic spectral correlators were evaluated by means of random matrix 
theory, and the universality of the result in that framework was established [||. The above expressions 
for the same quantities are more compact and convenient. We have explicitly checked in some special 
cases that the results of ref. || agree with those presented here. 

It is interesting to note that not only can we compute the microscopic spectral correlators directly 
from the corresponding finite-volume partition functions, we can also derive the universal double- 
microscopic limits of the orthogonal polynomials from these partition functions. This despite of the 
fact that these orthogonal polynomials seem to have no clear interpretation in field theory language. 

To derive expressions for the orthogonal polynomials, we make use of a convenient representation of 
these polynomials in terms of matrix integrals (see, e.g., ref. 

p£'H\;mi,...,m Nl ) = / Y[[dXMXi)] det^Af 1 ! [T( A " A ( 28 ) 

z 2n (mi,...,m Nf ) J -°° i= i i= i 



HidXM^)] det^Af 1 , (29) 

-00 ■ 



Here, 



and the measure factor w(Xi) is that of eq. (y). The relation fl28|) is readily verified by noting that it 
yields the required orthogonality relation. It also follows from ( p8| ) that the normalization corresponds 
to monic polynomials, i.e., P n (X) = X n + . . .. 

One sees that for n, N — * 00, with t = 2n/N fixed, eq. (pg|) becomes a relation between the orthogonal 
polynomials and (matrix model) partition functions, now involving an odd number of fermion species 
Nf + l (the additional fermion having mass iX). Let us now consider the double- microscopic scaling 
limit, and take t = l. The relation (28) then gives 

p 2 J ..,m f ) = c, z{Nf){ ^ — ^ , (30) 

where we again have replaced the matrix model partition functions by those of the finite-volume field 
theories (at the cost of introducing one overall proportionality constant Ci). 

We have indicated the relation for the even polynomials only. This is actually all that is required, 
since we can construct the odd polynomials from the following relation |B|: 

D (iV»/, v p2nil(X;m 1 ,...,m N )-P^ f \X;m 1 ,...,m N ) 
P 2n i 1 {X;m 1 ,...,m Nf ) = J — J — , (31) 
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where the even polynomials P^^A; mi, ■ ■ ■ , mN f ) are those of eq. (|2^), but in a different normaliza- 
tion: P 2 n (0) m l> • • • ' m Nj) = 1- 

The distinction between odd and even polynomials in the manner shown above actually has a consistent 
interpretation in terms of the finite-volume partition functions from field theory. Verbaarschot and 
Zahed have given these finite- volume partition functions in the form of integrals: 

Z (N J +1 \M) = JdU cosh[Tr{MUTU ] )] (32) 

for iV even, and 

Z( N f +1 \M) = JdU smh[Tr(MUTU^} (33) 

for N odd (where N denotes the three- volume) . Here M is the mass matrix already rescaled by NT,, 
whose diagonal form is 

M = diag(/i,/^i,/x 2 , •••,MAr / /2,-Mi,-A t 2,---,-MAr / / 2 ) , (34) 

and r = diag(ljv / /2+ii ~~ Ijvv/a)- We can choose to use just the even-iV partition functions p2] ) to get 
the even polynomials (|30[), and then derive the odd polynomials from the relation ([H]). Alternatively, 
we can use the odd-iV partition function directly (the analogue of the formula ( p8[ ) for odd polyno- 
mials). The result will be the same, as follows by expanding the partition function (32) for N + 2 
in terms of the partition function for N and then comparing with eq. ([H]). In fact, just from the 
relations (|3l| ) and (|30| ) it follows that if the even-iV partition function is given by (32), then the odd-iV 
partition function must necessarily be given by (|33|). 

Both of the above partition functions are given in terms of integrals of the SU (Nf + l)-mvariant Haar 
measure dU. These integrals can actually be evaluated rather easily by again making use of the 
Harish-Chandra-Itsykson-Zuber integration formula [13|. The case of one flavor of course stands on 
a special footing, and there one has, trivially, 

Z {1) (fi) = cosh(//) (35) 

for iV even, and 

£ (1) (/i) = sinh(/i) (36) 

for N odd. Here \x is the rescaled mass: fj,=mT*N, and we have dropped all irrelevant (/i-independent) 
overall constants. For higher values of Nf + 1 we find 

Z^Xin-m}) = ^i^[detD({ / ,; / ia) + (-l)^// 2 detD({- W -«})] , (37) 

where A(«M) is the the Vandermonde determinant of A4 eq. ( ^4|) and the additional sign (— \) N fl 2 
originates from A(— A4). The matrix D is of size (Nf + 1) x (Nf + 1) and is given by 

Dij = fjp- 1 ^ for 1 < j < — f - + 1 

N f Af, 

Dij = {-^y~ — - 2 e'^ for -J- + 2<j<N f + l 

i-i ,, N f N f 

D^ = ^-iV*- 1 for 2<i<^- + l ; 1 < j < -± + 1 

■ N f n Nf Nf 

Dij = r-V*- 1 for2<t<-^ + l; -f + 2 < j < N f + 1 



8 



4-i Nf Nf 

Dij = Mi-ie^- 1 for -J- + 2<i<N f + l ; 1 < j < -± + 1 

N f AT. AT, 

Dij = (-^xy-—- 2 e-^ for -J. + 2<i<N f + l ; -J- + 2 < j < N f + 1 . (38) 

We are now ready to state the results. For the polynomials we get 



H/ /7 (C 2 + ^) det f A M) A -M) \ 

V A ({-M) A ({M) J 

where in matrix D eq. (|3£|) has been replaced by iQ and A is the matrix defined in eq. (0). The 
normalization constant C\ remains undetermined, but can of course in any case be chosen at will. 
In H universal expressions for the same polynomials have been derived from random matrix theory. 
In contrast to the situation for the correlation functions eqs. ( p4[ ) and ( p6| ) the expression for the 
polynomials obtained here are less compact that in Q. We have checked explicitly in some special 
cases that both results agree. 

To illustrate the simplicity of the derivation presented here, let us look at the special case of quenched 
fermions Nf = 0. From eqs. (35) and plf ) we can immediately read off the well-known asymptotic 



Cl (-l) N f/ 2 det D(R; m}) + (-1) N+N f/ 2 det D({-*C; -Mi}) 



behavior of the orthogonal polynomials in the microscopic limit: 

Pf£(0 = C lC os(C) , P 2 ( ? +1 (C) = C lS in(C) . (40) 

Although we needed to evaluate the partition function for an odd number of fermions Nf + 1 to 
obtain the polynomials for an even number Nf, we have otherwise refrained from providing all the 
corresponding odd-iVj results for the double-microscopic spectral correlators, as well as the formulas 
for the orthogonal polynomials with an odd number of fermions. In fact, the physical interpretation 
of the odd case is not as straightforward as that of the even case. For example, in the massless 
limit the spectral density as normally defined is only positive definite for even TV" |2|. Moreover, the 
orthogonal polynomial technique does not directly apply, due to the measure being odd under parity 
(in eigenvalue space) for odd N. Indeed, even the formula (^) for the orthogonal polynomials will 
not be valid in that case, due to the non-existence of normalizable orthogonal polynomials. However, 
both partition functions (|32| ) and (|33"|), are completely well-defined, and can of course be used for 
convenience to compute the orthogonal polynomials for even Nf, as we have just shown. 



3 The Chiral Unitary Ensemble 



For the Dirac operator spectrum, this case corresponds to the gauge group SU(N C ), N c > 3 with Nf 
fermions in the fundamental representation. The matrix model partition function reads, in the sector 
of topological charge v (for convenience we shall consider u>0 throughout), Q 



(mi, 



,m N 



N f 

I dWY[ det (iM + m f ) exp 



-^trV(M 2 ] 



(41) 



where 



M 



W f 
W 



(42) 
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where W is a rectangular complex matrix of size N x (N+v), which is integrated over with the Haar 
measure. The space-time volume V of the gauge theory is, in the large- N (and large-T^) limit identified 
with 2N. 

Introducing the eigenvalues Aj of the hermitian matrix W^W, Z v can be written 

Nf N / N f \ 

zl Nf \m u ...,m Nf ) = nK)/°°n Yl(Xi + m 2 f )e- NV ^ IdetyA}" 1 ! 2 . (43) 

f=i Jo i=i v /=i / 

We have ignored all unimportant factors that arise from the angular integrations. Since the partition 
function for Nf fermions in the sector of topological charge v is related in a simple way to the partition 
function of the same Nf fermions plus v additional fermions of zero mass (in the sector of zero topo- 
logical charge), one can in principle restrict attention to the z^ = sector. It is nevertheless worthwhile 
to point out that the whole analysis which leads to a relation between the double-microscopic spectral 
correlators and the finite-volume partition functions carries over to the case of v 7^ 0. This will lead 
us to very compact expressions for these double-microscopic spectral correlators. 

The necessary generalization of the previous analysis to the present case of a chiral unitary ensem- 
ble with measure (|43|) is straightforward. We are here interested in the spectral correlations of M- 
eigenvalues z% rather than those of W-eigenvalues \ = z 2 . Because the whole procedure is identical to 
that of the previous section, and because the case j/ = already has been worked out in detail 0, we 
shall be brief. The two-point correlator, the kernel, is 

K N (z,z ;mi,... ,m Nf ) = 

Nf N -i 
e-f (n* 2 )+^' 2 ))(^+2 H ^( z 2 + m})(z' 2 + m 2 f) £ P(z 2 )P(z' 2 ) , (44) 

/ i=0 

where Pi{z 2 ) are the orthonormal polynomials associated with the above matrix model. The kernel 
can now be expressed as a normalized random matrix integral: 



K N J (z,z;m 1 ,...,m Nf) 

Z v (mi, . . . ,mjv 



N f „™iV-l 



X I1K) / II UA 4 Ar(A i -z 2 )(A i -z' 2 )n(A 4 + m ? )e 

f J ° 8=1 V f 



2 A -NV(Xi) 



f 

Thus, in the large- N limit we have 



detyA} -1 



(45) 



K^ f ' u \z,z';m u ...,m Nf ) = e~T ^ ^ (-If JIz 1 \[ ^f(z 2 + m})(z' 2 + m 2 ) 

f 

Z v (mi, . . . ,mjsi f ,iz,iz ) 



zj> Nf \mi, . . .,m Nf ) 



(46) 



By means of the usual factorization property, all higher n-point spectral correlation functions are then 
also explicitly expressed in terms of the two matrix model partition functions Zp S ^ and Zy 1 ^. The 
spectral density corresponds to the two additional (imaginary) masses being equal, as in eq. (0). 
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We now turn to the double-microscopic limit in which C, = zN2np(0) and \x% = rriiN2ivp{Q) are kept 
fixed as N^oo. The prefactor exp[— (N/2)(V(z 2 ) + V(z' 2 ))] again becomes replaced by unity, and by 
identifying £ = 2irp(0), we can now compare with the field theory finite- volume partition functions. 
This gives us the master formula 



/ 



7 y 



(N f +2) 



(pi,...,p Nf ,i(,i(') 

Zv \p>l, ■ ■ -,PN f ) 



(47) 



where the partition functions are those of the finite- volume field theories. Similarly, for the double- 
microscopic spectral density, 



Ps {(;pi,...,PN f ) 



C 2 |CII1(C 2 + /^ 



Zv (Pi, 



Zv(jJ>l, ■ ■ -,fJ>Nf) 



(48) 



All double-microscopic n-point correlation functions are again given by the factorization formula (18). 



A simple example which illustrates how powerful the above relations can be is that of the quenched 
case Nf = 0. All we need is the finite-volume QCD partition function for two massive fermions of 
degenerate masses ipf(N'E). This was evaluated analytically already by Leutwyler and Smilga 0] and 
found to be, in their normalization, 



Z^ 2 \ip,ip) = l v (ipf 



iu+l 



(ifi)I„-i(i(i) 



(49) 



where I n (x) is the nth modified Bessel function. The corresponding denominator in eq. (|47|) is again 
an irrelevant constant which we can set to unity. This gives 



P f v \0 = c 2 (-iy\(\ U(o 2 - j v+1 (0Ju-i(d 



The matching condition (1£) yields 



c 2 = -(-ir 



and hence 



pi M (C) = dCI [j,(C) 2 -J, +1 (C)^-i(0 



which is the known result [ffl] . Furthermore, by the previous considerations (cf. eq. 
that the general case of Nf massless fermions simply is equivalent to a shift v - 
recover the general massless result M without any effort: 



1 



|C| JN f +v{Cf - JN f +v+l(0 J N f +v-l(0 



(50) 



(51) 

3()) we also know 
v + Nf. We thus 



(52) 



We now need the general analytical expression for the finite-volume partition function for this case, 
with Nf fermions of arbitrary masses [jl5| (see also ref. p^| ). It can conveniently be written @ 



a^W-.w,) = d -^0f- (53) 



where the NfxNf matrix A({pi}) is 

Aj = p^Iu+j-iiPi) , (54) 
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and A(// 2 ) again indicates the Vandermonde determinant of the fj,f. 

For the numerator of eq. (p7|) we need the (Nf + 2) x (Nf + 2) matrix A with two of the entries being 
imaginary. This means that 

Aj = i u (-(i) j - 1 Jv+ j -i((i) for i = l,2, (55) 



and otherwise (for i>3) as in eq. (|54|). For convenience we pull out a factor of (-1) from every second 
column of the matrix .A, and also the factor of i v from the first two rows. This yields an overall factor 
of (— l^+I^// 2 ] where [x] denotes the integer part of x. Thus, 

Ko 1 (Ci,C2;M1> • • • jpNf) = C 2 , : — : , (56) 

(Ci 2 -C!)n /v /(Ci 2 + M^)(C! + M^) detA 

where the (Nf + 2) x (Nf+2) matrix B is defined by 

Bij = (-in- 2 ) i -' L Iu+j-i(lM-2) for 3<i<JV / + 2, (57) 
and the Nt x Nt matrix A is as in (|54]) . Using the Bessel relation 



^[x n J n+m (x)] = i n J n+m _i(3;)-rra n l J n+m (x) (58) 
we find the corresponding double-microscopic spectral density: 

where the (Nf + 2) x (iVj + 2) matrix S is defined by 

By = ( J - 2 Ju +3 -2(() (60) 



and = Bij for i^l. The general n-point correlators follow from eqs. (|18| ) and (|56|). Using the 
matching condition ( |l9[ ) gives 

C 2 = (-i)-W2] j 

and everything is now determined. 



For v = the results (56) and (p9|) agree with what has recently been obtained by a direct computation 
in random matrix theory |17| . While no explicit expressions were given for the case v ^ in ref. 
||, it could in principle be extracted from the general formulae for v = Q by setting v fermion masses 
equal to zero in a theory of Nf+v fermions. The result done in that way should of course agree with 
our explicit formula given above. Indeed, we have managed to prove by induction that the compact 
formulas (56) and ( |59| ) also follow from the v = results given in ref. Q. 

As in the previous section for the unitary ensemble, we now show that also in the chiral case the 
universal limit of the orthogonal polynomials can be obtained directly from the finite- volume partition 
functions alone. The extension of the formula ( |30|) to this chiral case is straightforward. There is 
now no parity "quantum number" for the polynomials, and the universal double-microscopic limit 



12 



is therefore unique, independent of whether the polynomials are of odd or even order. The explicit 
formula reads 



Pr f ' 1/ \X;m 1 ,...,m Nf ) = X1/ V / J[[dXMK)} det^H ]^[(A — Aj) , (61) 

Zvjiimi, . . . ,m Nf ) ->o i= i i= i 



where now the weight function w(X) is given by 

Nj 



w {\) = X"H[(X + m})] e-^W , (62) 



and 

N 



(at \ * r°° n 2 

zin f) = IIK) / Ui dx M^)\ <i«''„a; 1 , (63) 

Also here one easily sees that the normalization is such that the polynomials f ' u \\; mi, . . . ,mjsr f ) 



are momc. 



In the limit n — > 00 and N —* 00 with i = n/A fixed, the relation ( |6l| ) determines the orthogonal 
polynomials in terms of the matrix model partition functions. By going to the double-microscopic 
scaling regime with t = 1 where we can make the identification with the field theory partition functions, 
this gives us the relation 

n (C ;ni,...,VN f ) = C 3 {-1) «) ^ , (64) 

Z v (pi, . . .,/lNf) 

where the normalization constant C3 is still undetermined and we have passed to scaled M-eigenvalues. 
This relation is just as in the ordinary unitary case. For the numerator we need the (Aj + 1) x ( Ay + 1) 
matrix A eq. ( |54[ ) with one imaginary entry, 

M 3 = i v {-Q j - x J v+j -i{Q (65) 

and otherwise as in eq. (|54|). In order to fix the constant C3 and to compare with || we choose the 
normalization P^ f ' u \o; pL\ 1 . . . , /tijv>) = l. We then obtain 

{N s ,u) (t2 v _ uT[ M/ detP 

^ (C^i,..,^) - C ll C 2 + /x 2 deU ( 66 ) 



where 



V Xj = {-Cy- l J v+j ^{0 for l<j<A / + l , 

V ij = ri-l^+j-l&i-l) for 2 < i < Nf + 1 ; 1 < j < A/ + 1 , (67) 



and 



= rflv+jbH) for l<i,j<N f . (68) 

For 1/ = the above expression matches to the result for the polynomials in || , which was derived from 
random matrix theory. Here, eq. (^6|) directly gives the results for v zero modes, which is equivalent 
to set v fermion masses to in ||. We have proven by induction that the results of ref. || lead to 
precisely the same formula for arbitrary v as shown above. 
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4 The Symplectic and Orthogonal Ensembles 



As follows from the general classification of universality classes [jl|], gauge group SU(2) and Nf fermions 
in the fundamental representation correspond in matrix model language to the orthogonal ensemble, 
while gauge group SU(N C ) with Nf fermions in the adjoint representation of the gauge group corre- 
spond to the symplectic ensemble. 



The symplectic and orthogonal matrix ensembles are somewhat more complicated from an analytical 
point of view due to the non-existence of simple orthogonal-polynomial methods for those cases. The 
closest one apparently can get is based on the so-called quaternion method, which can be phrased 
in terms of skew-orthogonal (as opposed to truly orthogonal) polynomials (T^j. In this chapter we 
shall consider what may be the closest symplectic and orthogonal ensemble analogues of the relations 
derived above for the unitary and chiral unitary ensembles. 



We begin with the chiral symplectic matrix ensemble, as this is the case for which we most easily can 
derive useful relations that connect finite-volume partition functions to the associated Dirac spectra. 
In the language of random matrix theory, the partition function for the chiral symplectic ensemble is 
as in eq. (|4l|), except that now the integration is over matrices W whose elements are quaternion real 
Q. In terms of the eigenvalues Aj of the hermitian matrix W^W, Z u can now be written (we follow 
the conventional normalization, where the symplectic matrix model potential is rescaled by a factor 
of 2 compared with the chiral unitary case): 



Z u (mi, . . .,m Nf ) 



IT 



rrij 



I n h A 



2v+l 



m)) 



-2NV(Xi) 



dety-A) 



(69) 



Our goal is now to find the closest analogues of the master formulas (|lq) and (47). We shall make 



good use of some general relations that have been derived by Mahoux and Mehta [12]. Throughout 
this section we will use their notation here. 



The problem we encounter is that the quantity that most closely corresponds to the kernel of the now 
skew-orthogonal polynomials now is a quaternion /^(Aj, Aj), which can be represented by a 2x2 matrix. 
The correlation functions of eigenvalues are then given by quaternion determinants det[/4(Aj, Xj)] m 
of the kernel f±(\i, Aj). We have not been able to express this kernel itself in terms of matrix model 
(and thus also finite-volume field theory) partition functions, but only the determinants of this kernel, 
which are real valued functions. This will directly give us the expressions for the correlators, where 
we display the eigenvalue density and the density-density correlator as examples. Going back to 



eigenvalues Z{ of the Dirac operator rather than Aj 
following way: 



the spectral density can be obtained in the 



Pi {z;mi, 



,m Nf 




+ m 2 f ) 



-2NV(z 2 )_ 



in 



V 



Z^\m u ...,m Nf ) 



4^+3 



N f 

IK 

/=1 



z'a + m 



V 



-2NV(z*)i z 2 



v2|4 



2i-2 



+ m)) 



-2NV(z 2 ) 



zi N ^\ mi , 



,m Nf ,{iz}) 



zi Nf \mi, . . .,m Nf ) 



(70) 
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In the first step we have made use of Theorem 1.2 of ref. |12J in the form 

dz p+ idet[f 4 (zi,Zj)}p + i = (N - p)det[f 4 (zi,z j )]p , (71) 



using their explicit expression for the partition function. We have slightly generalized the measure of 
ref. [^] here to include the massive fermions and zeromodes. In the second step we have replaced the 
integral over N—l eigenvalues by the matrix model partition function with 4 equal massive flavors of 
imaginary mass iz, ignoring their difference in the large- N limit. 

When replacing the matrix model partition function in the scaling limit by its mesoscopic field theory 
counterpart, as was done in the previous sections, we obtain the following relation: 

Ps (C;Mi»- ■ ■ .Miv>) = C 11 (C + /*/) Tna 1 • ( 72 ) 

f =i Z u (pi, . . . , fi Nf ) 

Unfortunately the needed finite- volume partition function of 4 or more fermions in the adjoint rep- 
resentation is not known at present and thus the eigenvalue density cannot be further evaluated yet. 
But since from matrix model calculations the eigenvalue density is known to be expressible in terms of 
integrals of Bessel functions (see also [19]) a relatively simple expression for the field theory par- 
titions functions should exist. (Simple analytical formulas are indeed known at present up to Nf = 2 
§§)• 

As a second example we derive in a similar way an expression for the density-density correlator: 
p^\z,z';{rn f }) = j^TI) det[U(z,z% 



\z 2 - z' 2 \ 4 (zz') 4 »+ 3 TT (z 2 + m 2 )(z> 2 + m 2 ) e -2*(V(* 2 W 2 )) 

f=i zi N f\{m f }) 

l II d Zi z?+ 3 I] + m}) e~ 2NV ^\z 2 - z 2 \\z' 2 - z 2 



' 2 l 4 lWV TT 4- m 2 .W 2 4- ml\ p -2N(V(z 2 )+V(z' 2 )) 



{-lf u \z 2 - z' 2 \\zz'f H {z 2 + m)){z' 2 + m)) e 
/=l 

ff /+4+4) (mi, • • • , rriN f , {iz}, {iz'}) 
zi Nf \mi, . . .,m Nf ) 



(73) 



and thus in the double- microscopic limit, 



p,s (C,C • • • ,m f ) = w|C -C I (CO 11 CC + M/XC +Pj) jwa 

f=l Zv (Ml,---, P-Nf) 

(74) 

In the unitary ensembles the knowledge of this two-point correlator is equivalent to knowing the kernel 
as well, since its connected part is the square of the kernel p C on(z,z') = —K(z,z') 2 . However, in the 
symplectic case the kernel f^Zi^Zj) is not a symmetric function and does not factorize. All higher 
correlation functions can actually be expressed by partition functions in an analogous way as above. 
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Finally we turn to the case of the orthogonal ensemble, for which the relevant partition function, when 
expressed in terms of eigenvalue integrals, reads 



Z v (mi,...,m Nf 



N f ~, N I , Nt 

poo y I v _x ' 

IIK)/ U[ dX i X I 3 U^i + rn})e 

f=i Jo i=i V f=i 



-fV(A 4 ) 



i-1 



detijXj 



(75) 



Since the orthogonal case can only be treated on the same footing as the symplectic case using 

he same proc 

det[/i(z,z)]i 



quaternions 12], we will follow the same procedure as above: 



N 



N 



f 



z v n (z 2 + mj) e 
/=i 



2\ „-fV(« a ). 



1/ 



zi Nf \mx, . . . ,m Nf ) 



N-l 



N f 



n I ^ z i n 



( z 2 + m 2) e -fV(^)| r 2 



i '\z — z. 



i=l 



2i-2 



(76) 



and 



i 



iV(iV- 1) 



det[/ 1 (z,z')] 2 



N f 



iff Urn)) 



I 2 /2 

Z — 2 



/=i 



oo JV-2 N } 



det^zf- 2 



(77) 



It is straightforward to go to the double-microscopic limit of these expressions. However, in this case 
the absolute value inside the eigenvalue integrals eqs. ( |76| ) and ( |77| ) prevents us from immediately 
identifying them with the matrix model partition function with additional masses in any simple way. 
But since analogously to the symplectic ensemble an expression for the microscopic density in terms of 
integrals of Bessel functions is known [jlj , similar relations to finite- volume partition functions should 
exist. 



In this context it is particularly interesting to consider those relations between the kernels of the 
chiral symplectic/orthogonal ensembles and the chiral unitary ensemble which very recently have 
been derived by §ener and Verbaarschot pOfl . According to our present viewpoint, these relations not 
only extend the universality proof of ref. || to these ensembles, but also provide completely surprising 
and non-trivial identities among finite-volume partition functions for different effective field theories. 
An understanding of these new relations between effective partition functions would be very desirable 
at this point. 



5 Conclusions 



In this paper we have systematically explored the relationship between spectral correlators of the 
Dirac operator in the double-microscopic scaling region, and the corresponding finite- volume partition 
functions. Based on relations that can be proven in random matrix theory, we have shown how to 
extract the universal properties from these finite- volume partition functions alone. The most powerful 
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results are the two master formulas for the unitary and chiral unitary ensembles. These allow for a 
complete determination of all double-microscopic spectral correlators for those cases in terms of finite- 
volume partition functions. One of the surprising conclusions is that one can also derive the universal 
limits of the orthogonal polynomials of random matrix theory from the associated field theory partition 
functions. 

The resulting relations are therefore not just of interest from the field theory point of view (where 
they indicate that the formulation in terms of large- N random matrix theory to some extent can 
be avoided), but also from the viewpoint of random matrix theory. Indeed, in section 2 we have 
illustrated this in another way by showing how the famous (bulk) sine kernel of the unitary ensemble 
can be derived neatly from a simple SU(2) "chiral lagrangian". Many other examples can surely be 
found. 

The cases most exhaustively solved are those of the unitary and chiral unitary matrix ensembles. 
Here all pertinent information for the orthogonal polynomials, their kernel and thus all correlation 
functions are derivable from the corresponding finite-volume partitions, suitably extended to include 
more fermionic species of imaginary masses. For the chiral symplectic random matrix ensemble we 
have chosen a different way. We have directly related the correlation functions to the corresponding 
field theory partition functions with 4-fold degenerate additional fermion species. These relations may 
turn out to provide the most easy analytical derivation of the involved quantities. For the chiral 
orthogonal case there is a highly suggestive relation, which, however, relates the spectral density to 
the partition function of a theory where the absolute value of the determinant of the Dirac operator 
enters. More work is required here, to either relate this partition function to the conventional one, 
or to construct a chiral lagrangian that would correspond to taking the absolute value of the Dirac 
determinant. 



It is also worthwhile to reconsider the notion of random matrix universality in this new light. Of 
course, there is no substitute for the complete mathematical proof ||. But we can gain substantial 
insight into the mechanism of universality by tracing the disappearance of the matrix model potential 
V(A) in the relevant expressions. If we, for example, return to eq. (||) of random matrix theory 
and the master formula (|i~6|), we see that in the double-microscopic scaling limit the prefactor of 
exp[— (N/2)(V(x 2 ) + V(x' 2 ))] simply becomes replaced by unity. This is the only place where the 
random matrix theory potential enters explicitly (and where it disappears in the scaling limit, leaving 
a universal result). Corresponding factors of the potentials disappear in the other analogous relations. 
One should not be misled by these simple observations to conclude that universality of the matrix 
model results can be understood in such simple terms alone. To some extent the notion of universality 
is simply built into the crucial identification between matrix model and field theory partition functions 
in the mesoscopic, or double-microscopic, scaling regime. Indeed, the disappearance of factors such 
as exp[— (N/2)(V(x 2 ) + V(x' 2 ))] in the appropriate scaling limit is not the sole mechanism behind the 
proven universality of random matrix theory results. An obvious counterexample is provided by the 
recent study of microscopic limits of random matrix theories for which the macroscopic spectral density 



p(0) at the origin precisely is vanishing [21|. Here exponential prefactors such as those discussed above 
do not approach unity in the microscopic scaling limits, and still universality can be proven fl2l| . It 
would be interesting to find the chiral lagrangian analogues of these multicritical cases (for which 
p(O) = O0), and explore the relations discussed here in that more general context. 

2 There are hence no Goldstone modes, and no obvious group manifold on which to base the effective Lagrangian of 
the lightest hadronic excitations. A new principle seems needed to derive the corresponding effective Lagrangian. 
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